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PREFACE 


This  final  report  describes  the  results  of  work  performed  to  develop  a  model  of  noise  levels 
within  a  multi-engine  turbo-prop  aircraft.  The  work  was  performed  from  June  2001  to 
September  2002.  The  Principal  Investigator  for  this  project  was  Brian  Scott  Davis  of  the  Air 
Force  Research  Laboratory’s  Aural  Displays  and  Bioacoustics  Branch  at  Wright-Patterson  Air 
Force  Base  OH.  The  work  was  performed  under  Program  Element  62202F,  Work  Unit 
71841609.  Robert  A.  Lee  was  the  project  manager. 


m 


This  page  intentionally  left  blank. 


iv 


TABLE  OF  CONTENTS 


Page 


Section  1  -  INTRODUCTION . 1 

Section  2  -  EXPERIMENTAL . 4 

Section  3  -  THE  ANALYTICAL  MODEL . 17 

BACKGROUND . 17 

MODEL  A  (Monochromatic  Point  Source) . 21 

MODEL  B  (Extended  Source  Function) . 26 

MODEL  C  (Uniformly  Rotating  Paddle) . 35 

Section  4  -  DISCUSSION  And  CONCLUSIONS . 43 

References . 46 


v 


LIST  OF  FIGURES 


Figure  Page 

2.1.  C- 130  Scale  Top  View . 4 

2.2.  C- 130  Scale  Frontal  View . 5 

2.3.  Sampled  Microphone  Fuselage  Locations . 8 

2.4.  One/Third  Octave  Band  Spectra  at  Microphone  Seven  - 

Analog  Synchrophaser . 12 

2.5.  One/Third  Octave  Band  Spectra  at  Microphone  Seven  - 

Digital  Synchrophaser . 12 

2.6.  Linear  Interpolation  of  the  Measured  Data  - 

Phase  Setting  (180°,  180°,  180°,  180°) . 13 

2.7.  Linear  Interpolation  of  the  Measured  Data  - 

Phase  Setting  (150°,  180°,  150°,  120°) . 14 

2.8.  Linear  Interpolation  of  the  Measured  Data  - 

Phase  Setting  (180°,  180°,  150°,  120°) . 15 

2.9.  Linear  Interpolation  of  the  Measured  Data  - 

Baseline  Analog  Controller  Configuration . 16 

3.1.  Predicted  Fuselage  Acoustic  Intensity  -  Model  A  - 

Phase  Angle  (180°,  180°,  180°,  180°) . 23 

3.2.  Predicted  Fuselage  Acoustic  Intensity  -  Model  A  - 

Phase  Angle  (150°,  180°,  150°,  120°) . 24 

3.3.  Predicted  Fuselage  Acoustic  Intensity  -  Model  A  - 

Phase  Angle  (180°,  180°,  150°,  120°) . 25 

3.4.  Pulsating  Disk  Propeller  Assembly  Representation . 26 

3.5.  Predicted  Fuselage  Acoustic  Intensity  -  Model  B  - 

Phase  Angle  (180°,  180°,  180°,  180°) . 32 

vi 


LIST  OF  FIGURES  (Continued) 


3.6.  Predicted  Fuselage  Acoustic  Intensity  -  Model  B  - 

Phase  Angle  (150°,  180°,  150°,  120°) . 33 

3 .7.  Predicted  Fuselage  Acoustic  Intensity  -  Model  B  - 

Phase  Angle  (180°,  180°,  150°,  120°) . 34 

3.8.  Predicted  Fuselage  Acoustic  Intensity  -  Model  C  - 

Phase  Angle  (180°,  180°,  180°,  180°) . 40 

3.9.  Predicted  Fuselage  Acoustic  Intensity  -  Model  C  - 

Phase  Angle  (150°,  180°,  150°,  120°) . 41 

3.10.  Predicted  Fuselage  Acoustic  Intensity  -  Model  C  - 

Phase  Angle  (180°,  180°,  150°,  120°) . 42 


vii 


LIST  OF  TABLES 


Table  Page 

2.1.  Measured  Phase  Angle  Combinations . 6 

2.2.  X-Y  Coordinates  of  Microphone  Locations . 9 


viii 


Section  1  -  INTRODUCTION 


The  C-130  is  a  four-engine  turbo-prop  air  vehicle  principally  designed  to  ferry  military 
payloads  in  excess  of  45,000  pounds.1  High  levels  of  noise  generated  by  powerful  aircraft 
engines  is  one  of  the  by-products  of  being  able  to  airlift  heavy  military  cargo.  For 
example,  C-130  flight  crew  personnel  typically  experience  sound  pressure  levels  (SPL) 
between  100  and  115  dB,  15-25  dB  above  the  OSHA  standard  of  90  dB  per  eight-hour 
workday.2  Intense  noise  exposures,  unfortunately,  are  not  unique  experiences  and  often 
result  in  hearing  loss  for  many  military  personnel.  According  to  data  obtained  from  the 
Department  of  Veterans  Affairs,  the  total  cost  of  compensation  in  2001  was 
$361,362,000  for  all  veterans  having  hearing  loss  identified  as  a  primary  disability.3  This 
permanent  impact  on  human  hearing  and  monetary  cost  could  potentially  be  mitigated  by 
more  widespread  application  of  active  noise  cancellation  techniques. 

Active  noise  cancellation  involves  creating  an  equal  amplitude  but  opposite  phase 
waveform  to  cancel  an  acoustic  source.  Area  cancellation  techniques  within  an  aircraft 
typically  involve  coating  the  walls  of  the  fuselage  with  lightweight  acoustic  damping 
material.  Although  this  approach  works  well  for  jet  aircraft  where  a  significant  amount 
of  the  noise  generated  by  the  engines  is  above  500  Hz,  it  has  little  effect  on  turbo-prop 
aircraft  where  the  spectra  is  dominated  by  low  frequency  tones  associated  with  the  blade 
passage  frequency  (BPF).4,5  Control  of  these  longer  wavelengths  is  generally  expensive 
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and  difficult  since  passive  noise  control  measures  require  large  muffler  assemblies  or 
heavy  enclosures.4,6  Active  noise  technology  provides  the  solution.  The  commercial 
sector  has  already  implemented  this  technology  to  minimize  the  noise  generated  from  air 
handling  system  exhaust  ducts,  automobile  engine  vibration,  and  low  frequency  tonal 
noise  from  mid-sized  propeller  aircraft.6 

Successful  reduction  of  the  acoustic  levels  within  the  fuselage  of  turbo-prop  aircraft  has 
been  previously  demonstrated  using  two  distinct  methods.  Method  One  utilizes  dedicated 
sensors  positioned  in  the  vicinity  of  each  passenger’s  head  to  sample  the  energy 
generated  by  the  propellers  at  the  Blade  Passage  Frequency  (BPF)  and  its  harmonics. 

This  data  is  processed  in  real  time  to  produce  an  error  signal  input  to  an 
amplifier/loudspeaker  system  capable  of  producing  a  secondary  field  which  is  180 
degrees  out  of  phase  with  the  input.  Although  effective,  installation  of  multiple  amplifier 
and  speaker  arrays  capable  of  producing  a  signal  of  the  required  amplitude  (100-1 15  dB) 
within  the  fuselage  is  cost  and  space  prohibitive.5, 7’ 8 

In  an  alternative  approach,  cancellation  is  accomplished  at  the  source  prior  to  entry  into 
the  fuselage  by  enhanced  control  of  the  relative  phase  between  propellers  (Method 
Two).9  Microphones  are  again  utilized  to  develop  an  error  signal  that  functions  as  an 
input  to  an  advanced  propeller  phase  controller,  known  as  a  synchrophaser.  The 
synchrophaser  couples  the  error  signal  with  a  once  per  revolution  signal  of  the  propeller 
engine  to  issue  speed  and  phase  control  commands  to  the  propeller  governor.  Propeller 
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speed  and  phase  are  maintained  through  this  method.  Phase  angles  are  either  commanded 
manually  or  automatically  via  an  adaptive  search  routine. 

The  purpose  of  this  thesis  is  to  develop  an  active  noise  cancellation  analytical  model 
(following  Method  Two)  that  can  accurately  predict  the  acoustic  intensity  field  within  the 
fuselage  of  a  turbo-prop  aircraft  with  respect  to  the  relative  propeller  phase  angle 
between  the  four  engines.  Theoretically,  the  physical  model  would  enable  the 
development  of  engines  and  control  systems  specifically  designed  to  minimize  interior 
noise  independent  of  airframe  design. 
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Section  2  -  EXPERIMENTAL 

The  test  vehicle  utilized  for  these  acoustic  tests  was  an  operational  United  States  Air 
Force  MC-130E.  The  C-130  is  an  all-weather  cargo  transport  specifically  tasked  to 
perform  infiltration  and  exfiltration  operations  in  hostile  environments.  Testing 
conducted  in  January  2002  spanned  three  weeks  at  Eglin  Air  Force  Base,  Florid  During 
these  three  weeks,  25  flight  hours  of  data  was  recorded.  Schematic  drawings  of  the 
aircraft  dimensions  are  provided  in  Figure  2.1  and  Figure  2.2. 10 


Figure  2.2  C-130  Scale  Frontal  View.  Engines  are  numbered  1-4  from  right  to  left. 

Four  Allison  turbo-prop  engines  with  an  effective  shaft  horsepower  of  4,300  power  the 
aircraft.  Each  of  these  turbojet  engines  drives  a  four  bladed  Hamilton  Standard  propeller 
assembly  that  maintains  a  constant  speed  of  1020  revolutions  per  minute  (RPM).1’ 10  This 
corresponds  to  a  primary  BPF  of  68  Hz  and  an  associated  wavelength  of  5  meters. 
Propeller  direction  is  counter  clock  wise  with  respect  to  an  observer  situated  forward  of 
the  engines.  Engines  were  identified  1-4,  right  to  left  (left  to  right  with  respect  to  the 
pilot  position).  Thus  in  Figure  2.1  and  Figure  2.2,  engine  2  is  the  inboard  engine  to  the 
right  of  the  fuselage.  Inter-engine  distances  of  5  meters  for  the  inboard  and  10  meters  for 
the  outboard  were  calculated  based  upon  scale  diagrams  Figure  2.1  and  Figure  2.2.  All 
values  are  measured  with  respect  to  a  coordinate  system  situated  at  the  intersection  of  the 
line  formed  by  the  propeller  hubs  and  the  plane  that  bisects  the  fuselage  from  nose  to  tail. 

Because  the  C-130  is  an  aircraft  with  synchrophased  engines,  an  analog  controller  has 
routinely  been  utilized  to  maintain  a  fixed  phase  relationship  between  the  engines.  Either 
of  the  two  inboard  engines  can  be  identified  as  the  master  reference.  Phase  relationships 
were  maintained  with  respect  to  the  master  and  changes  were  accomplished  via  small 
perturbations  in  propeller  RPM.11  Since  the  phase  angles  were  invariant,  the  analog 
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system  did  not  compensate  for  variations  in  the  sound  field  as  a  function  of  aircraft 
velocity.  This  resulted  in  poor  management  of  the  interior  noise  profile. 
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Table  2.1  Measured  phase  angle  combinations.  Engine  #2  was  selected  as  master  for  all 
measurements.  The  phase  angle  for  this  engine  was  arbitrarily  set  at  180°  because  the  digital 
controller  could  not  stabilize  at  zero  degrees.  Highlighted  combinations  identify  particular 
combinations  for  comparison  to  model  results  from  the  next  section. 


For  the  acoustic  test  here,  a  form-fit  function  replacement  was  used  (along  with  the 
original  analog  controller),  based  on  an  active  synchrophaser  prototype  developed  by 
British  Aerospace  Electronic  (B.A.E.)  Systems.  In  this  design,  the  user  maintains  the 
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ability  to  select  a  master  engine  and  gains  the  ability  to  dial  in  a  particular  set  of  phase 
angles.  These  capabilities  facilitate  better  control  of  the  noise  signature  of  the  aircraft. 

A  two-tiered  data  collection  approach  was  adopted  to  meet  the  funding,  manpower,  and 
aircraft  availability  constraints.  The  first  tier  consisted  of  acoustic  measurements  of  both 
the  internal  and  external  sound  fields  at  27  sets  of  phase  angles.  The  phase  angles 
measured  in  tier  one  are  identified  in  Table  2.1.  Engine  2  was  designated  as  the  master 
engine.  All  other  engines  were  varied  with  respect  to  the  fixed  orientation  of  this  engine. 

The  primary  objective  of  this  particular  test  was  to  localize  the  phase  angle  set  that 
produced  a  significant  reduction  in  the  external  sound  field;  these  data  sets  were  to  be 
analyzed  in  real  time  such  that  a  crude  estimate  could  be  made  as  to  where  the  global 
minimal  set  might  reside.  A  second  set  of  phase  angles  (tier  two  data)  was  then  to  be 
performed  to  attempt  to  determine  the  minimal  set.  However,  analysis  of  the  external 
data  in  the  field  proved  to  be  inconclusive  and  the  tests  were  terminated  after  the  tier  one 
data  was  acquired  for  two  airspeeds,  1 15  and  240  knots.  These  airspeeds  represent  two 
typical  operational  flight  conditions.  The  analysis  of  the  internal  data  is  of  primary 
interest  in  this  thesis.  Further  analysis  of  the  external  data  is  the  subject  of  ongoing 
research  at  the  Air  Force  Research  Laboratory. 

The  system  utilized  to  collect  the  acoustic  information  consisted  of  sixteen  Acousti-Cell 
microphones,  a  power  supply,  and  a  32-channel  METRUM  digital  tape  recorder.  The 
microphones  were  secured  to  a  cargo  strap  connected  between  the  paratroop  anchor 
cables  that  extend  the  length  of  the  fuselage  cabin  and  positioned  slightly  above  head 
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height  as  to  not  interfere  with  normal  crew  operations.  In  addition  to  the  raw  pressure 
field  data  from  the  16  microphones,  channels  17-20  recorded  the  once-per-revolution 
signal  from  each  engine  to  verify  that  phase  control  was  maintained  during  measurement. 
Channel  21  was  reserved  for  an  IRIG-B  time  code  signal  to  synchronize  the  data  upon 
post-analysis.  Unfortunately,  the  32-channel  METRUM  recorder  failed  during  system 
tuning  and  a  TEAC  RX-816  16  channel  digital  tape  recorder  had  to  be  installed  in  its 
place.  With  only  16  channels  available,  five  microphone  channels  had  to  be  sacrificed. 
The  1 1  channels  sampled  are  identified  in  Figure  2.3. 


Figure  2.3  Sampled  Microphone  Fuselage  Locations.  FS  stands  for  fuselage  station  and  refers  to  the 
distance  measured  from  the  nose  of  the  aircraft  in  inches. 

These  microphone  locations  were  selected  in  order  to  prevent  the  creation  of  significant 
holes  in  the  data  set.  The  positions  of  the  cargo  bay  microphones  in  a  y-z  coordinate 
system  centered  on  the  axis  of  symmetry  (z-axis)  and  the  line  joining  the  propeller 
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centers  (y  -axis)  are  represented  in  Table  2.2.  Information  from  cockpit  microphones  14 
and  16  will  be  ignored  since  the  primary  concern  is  cargo  bay  noise  levels.  It  was 
assumed  that  microphones  8-11  resided  in  the  propeller  plane  of  rotation  since  the 
separation  was  much  less  than  one  wavelength.  Signals  were  recorded  once  per 
revolution  on  channels  1,  3, 6, 10  and  the  IRIG-B  signal  was  recorded  on  channel  16. 
The  PC  controller,  power  supply  and  other  recording  equipment  were  positioned  at  T1 
and  T2. 


Mic  Number 

y  coord(m) 

z  coord  (m) 

8 

-1.07 

0.00 

9 

-0.53 

0.00 

11 

1.04 

0.00 

2 

0.00 

-6.53 

4 

-1.07 

-2.34 

6 

0.53 

-2.34 

7 

1.07 

-2.34 

12 

-0.53 

2.54  i 

13 

0.53 

2.54 

Table  2.2  Y-Z  Coordinates  of  Cargo  Bay  Microphone  Locations. 

Each  microphone  was  calibrated  using  a  1 14  dB,  250  Hz  tone  before  each  test  flight  and 
the  respective  Engineering  Unit  (EU)  scale  factors  for  each  microphone  were  recorded  to 
facilitate  post-analysis.  All  microphone  operations  were  validated  against  a  gold  standard 
microphone  calibrated  by  the  National  Institute  of  Standards.  Laboratory  analysis  of  the 
recorded  data  was  conducted  using  two  HP35665A  Dynamic  Signal  Analyzers  (DSAs). 
The  DSAs  were  configured  to  process  the  raw  data  over  a  frequency  range  of  10  Hz  to  16 
kHz  referencing  the  EU  scale  factor  recorded  in  the  field  for  each  microphone.  Both 
analyzers  were  set  to  a  0.25  integration  time.  A  typical  one-third-octave  band  spectrum 
for  the  C-130  in  standard  analog  controller  configuration  (phase  angles  unknown)  is 
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shown  in  Figure  2.4  for  Microphone  seven  and  a  similar  plot  at  the  same  location  is 
shown  in  Figure  2.5  for  the  digital  control.  (Phase  Angles  150°,  180°,  150°,  180°). 

In  Figure  2.4,  the  contributions  of  the  BPF  at  108  dB  are  extremely  evident  as  the  63  Hz 
band  is  approximately  7  dB  above  the  background  jet  turbine  noise.  Information  from 
the  68  Hz  BPF  appears  in  the  63  Hz  band  under  one/third  octave  band  analysis.  Whereas 
in  Figure  2.5  with  the  digital  synchrophaser  engaged,  the  BPF  blends  in  with  the 
background  radiation  at  a  particular  phase  setting  of  150°,  180°,  150°,  180°.  There  are  no 
observable  higher  harmonics;  if  present  they  are  buried  within  the  background  radiation. 
The  variance  in  the  data  as  a  function  of  time  for  the  enhanced  controller  was  generally 
less  than  3  dB.  Higher  variances  occurred  on  Microphones  4  and  6  near  the  conclusion 
of  the  flight  test.  These  variances  were  attributed  to  systematic  microphone  failure.  The 
data  from  each  microphone  for  a  given  phase  condition  was  plotted  as  a  function  of 
position  and  linearly  interpolated  using  MathCAD  to  produce  contour  plots  of  the 
acoustic  sound  field  within  the  fuselage. 

Three  particular  phase  angle  sets  were  of  interest.  The  first  set  (180°,  180°,  180°,  180°) 
was  the  condition  that  should  produce  constructive  interference.  The  sound  field  for  this 
set  is  plotted  in  Figure  2.6.  The  second  is  a  combination  that  reflects  destructive 
interference.  Inspection  of  the  plots  indicated  that  a  higher  degree  of  destructive 
interference  occurred  when  the  inboard  engines  are  only  30  degrees  out  of  phase.  The 
phase  set  chosen  for  comparison  (150°,  180°,  150°,  and  120°)  is  plotted  in  Figure  2.7. 

The  last  set  chosen  for  comparison  was  the  phase  set  that  seemed  to  subjectively 
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produced  the  calmest  distribution  within  the  fuselage:  (180°,  180°,  150°,  120°).  See 
Figure  2.8.  The  baseline  analog  controller  is  illustrated  in  Figure  2.9.  Actual  aircraft 
flyover  photos  of  the  aircraft  at  each  phase  setting  provided  verification  of  phase  lock  and 
are  compliments  of  MSgt  Alan  Martin,  AFRL/HECB.  All  plotted  intensities  are  at  the  63 
Hz  1/3  octave  band  (dB)  and  are  plotted  against  a  1/90  scale  aircraft  schematic.1  Cockpit 
microphone  data  was  not  included  in  the  plots. 
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Figure  2.6  Linear  Interpolation  of  the  Measured  Data  -  Phase  Setting  (180°,  180°,  180°,  180°)  at 
63Hz,  One/Third  Octave  Band  (dB)1  -  240  Knots 
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90  100  110  120 

Figure  2.7  Linear  Interpolation  of  the  Measured  Data  -  Phase  Setting  (150°,  180°,  150°,  120°)  at 
63Hz,  One/Third  Octave  Band  (dB)  -  240  Knots 
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Figure  2.8  Linear  Interpolation  of  the  Measured  Data  -  Phase  Setting  (180°,  180°,  150°,  120°)  at 
63Hz,  One/Third  Octave  Band  (dB)  -  240  Knots 


90  100  110  120 

Figure  2.9  Linear  Interpolation  of  the  Measured  Data  -Baseline  Analog  Controller  Configuration  at 
63Hz,  One/Third  Octave  Band  (dB)  -  240  Knots 
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Section  3  -  THE  ANALYTICAL  MODEL 

BACKGROUND 


In  order  to  develop  an  analytical  model  for  the  acoustic  pressure  fluctuation,  A P,  at  a 
given  point,  Q,  inside  the  fuselage,  we  must  first  develop  a  mathematical  expression.  The 
most  obvious  starting  point  is  with  the  complete  lossless  inhomogeneous  acoustic  wave 
equation,12 


2  l  32/>  _  dG  -  a  (pujuj) 
c2  dt2  dt  dxfy  ’ 


where  c  is  the  wave  speed,  G  is  the  rate  at  which  mass  is  injected  into  a  volume,  F  is  the 


^2 

body  force  per  unit  volume,  and  is  the  contribution  from  viscous  forces. 


dxjdxj 


This  first  term  in  equation  (3.1)  above  can  be  eliminated  as  our  propeller  source  is  a  rigid 
body  and  is  not  ejecting  matter  into  the  fluid.  Additionally,  the  last  term  can  also  be 
eliminated  since  our  source  is  not  the  product  of  spatial  rates  of  change  within  the  fluid  as 
in  the  exhaust  of  a  jet  engine.  Thus  the  only  term  that  survives  is  the  middle  term,  which 
represents  body  forces  as  those  forces  produced  by  moving  through  a  fluid  without  a 
change  in  volume.  So  we  have  an  inhomogeneous  equation  of  the  following  form 

V2P-4^  =  V-F.  (3.2) 

c2  dt2 

This  equation  has  the  basic  structure. 


V2xP(r,  t)-49  =  -4 nfCr,  0 , 
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where  f(r,  t)  is  some  source  function  that  includes  a  spatial  and  time  dependence.  It  is 
often  customary  to  solve  an  equation  of  this  form  by  developing  a  Green’s  function  that 
expresses  how  the  source  disturbance  will  propagate.13 

To  develop  the  propagator  for  this  case,  we  first  assume  that  we  have  a  simple  system 
devoid  of  boundary  surfaces.  Under  these  conditions,  the  explicit  time  dependence  can 
be  eliminated  by  using  the  Fourier  integral  transforms, 

oo 

'P(r,  t)  =  —  |V(F,  (D)exp(-i(Ot)  d CO 
2nJ 

OO 

f(r,t)  =  -J-  f  /(?,  (0 ) exp(-MMf)  d co .  (3.4) 

2n  J 

— oo 

By  direct  substitution  of  these  integrals  into  our  basic  equation  for  the  acoustic  field  (3.3) 

we  arrive  at  the  equation, 

2 

V2^?,  c o)+^'¥(i ,  (O)  =  -An f{r,  a) ,  (3.5) 

c 

Equation  3.5  is  merely  the  nonhomogeneous  scalar  Helmholtz  equation,  where  G)/c  =  k, 
the  wave  vector  associated  with  the  frequency  to. 

The  Green’s  function  associated  with  equation  (3.5)  is 

(v2  +  k  2)Gk(r,r)  =  -4nS(7-r).  (3.6) 

If  boundary  surfaces  are  not  present,  the  surface  integral  in  the  Green’s  function  vanishes 
and  Gk (r,r)  depends  only  on  the  distance  between  the  field  point  and  the  source|r  -  ?\ . 

The  Laplacian  operator  requires  that  Gk  must  satisfy. 
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(3.7) 


±.^{RGk)+k^Gk=-AnS(R), 

where  R  represents  the  distance  from  the  source  point  to  the  observation  location,  Q. 
Hence,  we  have  a  general  solution  that  is  defined  for  all  points  in  space  except  where  R  = 
0. 

RGk  ( R )  =  A  exp (ikR)  +  Bexp(-ikR)  (3.8) 

In  the  limit  kRQ  l  and  the  equation  reduces  to  Poisson’s  equation.  Thus  the 

general  normalized  solution  to  the  time-independent  Green’s  function  is 

q  ^  _  Aexp(ikR)  +  Bzxp(-ikR)  ^9) 

R 

with  A  +  B  =  1.  Now  that  the  time-independent  equation  has  been  developed,  we  need 
to  derive  the  time-dependent  version.  (The  time  independent  case  will  be  explored  in 
Model  A.) 


The  corresponding  Green’s  function  for  the  time  dependent  case  must  satisfy 

V~L  4L)Gk{r,t-,?S)  =  aS{r-r')S{t-t')  (3.10) 

l  c  dt  ) 

The  time  Fourier  transform  of  Gk  (r,r,f\t')c an  be  shown  to  be  Gk  (R)exp(i(ot’) .  Therefore, 

oo 

Gk  (r,t;r'S)  -  —  f  Gk{R)t^p(icot^ty^{-iQ)t)dQ) , 

2n  J 


Gt(r,/;r',r')  =  G*vfr-— 1  +  G»_<yfr+— ), 

V  c )  \  c ) 

Gk(r,r,r,t')  =  Gk+{R,T)+Gk-{R,T), 


(3.11) 
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where  x-t-t' . 


Expressing  the  propagator  in  an  explicit  form, 


Gki(R,r)  = 


_  v 


r-r | 


r-AW 


JJ 


IM 


(3.12) 


We  neglect  the  advanced  Green’s  function  Gk~(R,x)  because  we  are  only  concerned  with 
the  particular  delta  function  argument  that  describes  the  acoustic  field  at  some 
observation  point  Q  at  time  t  caused  by  a  disturbance  at  some  distance  r  away  at  an 
earlier  time  t ' .  Gk+(R,x)  is  sometimes  called  the  causal  Green’s  function.14  Now  with  the 

propagator  completely  defined,  we  can  determine,  at  least  in  principle,  the  acoustic  field 
at  any  point  by  integrating  the  product  of  our  Green’s  function  and  the  source  function  , 

S,  over  all  space  and  time,  neglecting  surface  effects: 

AP(r,t)=  jjGk+ (r,f,r',t')S(r',t)  dr'dt' .  (3.13) 


At  this  point,  the  model  development  becomes  a  little  more  complicated  as  the 
wavelength  (A  =  5m)  is  on  the  same  order  as  |F  -7  j .  This  region  of  interest  is  generally 

referred  to  as  the  intermediate  or  induction  zone15  and  makes  evaluation  of  the  integral  in 
(3.13)  difficult.  The  neglect  of  any  surface  present  means  that  refraction  and  reflection 
effects  will  be  ignored  These  effects  are  considered  small  since  the  relevant  wavelength 
is  an  order  of  magnitude  larger  than  the  thickness  of  the  fuselage  skin. 


In  the  remainder  of  this  section,  three  models  will  be  developed  based  upon  the  foregoing 
discussion. 
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MODEL  A  -  Monochromatic  Point  Source 

In  our  first  attempt  at  modeling  the  system,  we  shall  take  a  simple  case  that  should  be 
relatively  easy  to  evaluate,  namely,  treating  the  engines  as  point  sources  having  a 
harmonic  time  dependence  exp  i(o»-^(  )  where  a)  is  the  BPF  (68Hz)  and  (J);  is  the  phase 
angle  of  the  Ith  engine  (i  =  1-4) .  In  this  case,  the  source  term  for  the  Ith  engine  can  be 
written. 


Sj(r,t)  =  StS(r  -  r')exp(-i(cot  -  $ ))  (3.14) 

where  the  constant  Sj  denotes  the  acoustic  strength  for  propeller  noise  associated  with 
engine  /  centered  at  r  From  equation  (3.13)  we  obtain  the  effect  of  a  point  source  at 
engine  i  on  the  acoustic  amplitude  at  a  point  7  and  at  a  time  t. 


AIJM- 


Sj  exp(i0>  )exp(ifc  |r  -  /j  j)  . 
- ,  ' - exp(-jfitf). 


(3.15) 


where  jfc  =  — . 

C 


The  combined  effect  of  the  four  engines  is,  by  superposition, 
ZUJ  (r.Q  S‘  'Ll  p. 

tf  Ml 


(3.16) 


The  sound  intensity  /(F)  at  the  observation  point  7 ,  averaged  over  one  cycle,  is 
proportional  to, 


d/_  n2  _  ^  StS  *  exp( j(*  |F  -  /;  j +«?,)-/(*  |F  -  Fy  ]  +  (pj )) 
i[r’t}  =k  \7-rA\r-rA 


(3.17) 


=isl2  Zz4tt+2S 

hH  <» 


COS 

’(* 

r-r, 

1  )~{k7-7j 

\7-7i]7-7j\ 

(3.18) 


assuming  the  sources  are  all  of  equal  strength,  S. 
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Equation  (3.18)  was  used  to  calculate  the  relative  acoustic  intensity  at  an  array  of  points 
within  the  fuselage  of  the  aircraft.  Since  the  experimental  data  is  recorded  in  decibels 
(dB),  the  calculated  values  were  scaled  according  to, 


I(dB)  =  scale_factor  + 


lOlog 


f— Y 

.Vef  JJ, 


(3.19) 


where  Iref  was  taken  to  be  20  micropasc  a  and  the  scale  factor  was  chosen  to  make  the 
I(dB)  fall  within  the  range  of  the  experimental  data  set  (90  to  120  dB).  With  the  values 
calculated  for  the  acoustic  intensity  the  data  was  interpolated  in  the  same  manner  as  the 
actual  data  sets.  Representative  cases  are  included  in  Figures  3. 1-3.4. 
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Figure  3.1.  Predicted  Fuselage  Acoustic  Intensity  -  Model  A  (Point  Source)  -  Phase  Angle  (180°,  180°, 
180°,  180°). 


Figure  3.2.  Predicted  Fuselage  Acoustic  Intensity  -  Model  A  (Point  Source)  -  Phase  Angie  (150°,  180°, 
150°,  120°) 
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Figure  3.3.10  Predicted  Fuselage  Acoustic  Intensity  •  Model  A  (Point  Source)  -  Phase  Angle  (180°, 
180°,  150°,  120°) 


MODEL  B  -  Extended  Source  Function 


Consider  now  that  each  propeller  assembly  can  be  visualized  as  a  pulsating  disk  of  radius 
“a”,  constrained  to  the  plane  of  rotation  of  the  propeller  assembly.  See  Figure  3.4. 


Figure  3.4.  Pulsating  Disk  Propeller  Assembly  Representation 

In  addition  to  being  an  extended  rather  than  a  point  source,  the  time  dependence,  fit'), 
although  assumed  periodic  with  a  period  equal  to  the  BPF,  will  not  be  restricted  to  being 
harmonic.  Thus, 

Sir, t)  =  Six', y^Siz1)  fiO,  (3.20) 

where  the  primed  indices  denote  the  source  coordinates  and  the  Dirac  delta  function  in 
3.25  guarantees  no  z '  contribution.  If  fit')  is  non-harmonic  there  will  be,  in  addition  to 
the  fundamental,  o\  =  BPF ,  higher  harmonics  contributing  to  the  acoustic  signal. 
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We  assume  that  S(x',y')  is  radially  symmetric  in  the  x'y'  plane  about  the  propeller  hub. 
Then  inserting  the  source  function  into  3.18  and  integrating  over  z'and  t' yields, 


(3.21) 


Expressing  the  time  contribution  as  a  Fourier  series. 


fit)  =  £c„  exp(-MW„/ ) exp(i<p) , 


(3.22) 


where  Cn  is  the  n*  Fourier  component,  con  =  no  and  (J)  is  a  phase  angle,  AP(r,t)  can  be 


written  as, 


A P(r,t)  =  exp (i<p)*rcn  exp (-io„t)S„(r) , 


(3.23) 


where 


a  2tr 

Sn  (r  )  =  J  J  j|iLLexp(i*n  | r  -  r  '|)r  'dr  'd<D 
0  0  1  r 


(3.24) 


j  ,  no) 
and  k„  =  — , 


If  we  denote  by  the  index  “j”  the  contribution  from  the  j1*1  engine,  then  the  total  pressure 
fluctuation  at  the  observation  point  r  at  a  time  t  is 


AP(r,0  =  ^AP;(r,0 


(3.25) 


Multiplying  this  expression  by  its  complex  conjugate  and  taking  the  time  average  over  a 
period  T  gives 


| AP(r,r)|  =  ~ 2  2 exP ) exp H<Pr)Yj Sn - on)dt 

*  .'I  ill  _ 1  ‘  .  - 


n= 1  n  =1 
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By  virtue  of  orthogonality  ,  only  the  terms  with  n=n'  survive.  Hence, 


exp  <k>j)Sf(r) 
J=l 


(3.26) 


Now  all  that  remains  is  to  evaluate  the  Snj(7) .  We  shall  assume  that  S(r')  is  a  constant,  So 


Referring  to  Figure  3.4,  the  distance  between  a  field  point  7  and  a  source  point  f  can  be 
expressed  as  a  series  of  Legendre  polynomials, 


r~r\  Vrj 


(3.27) 


where  y  is  the  angle  between  7  and  7’ .  The  series  converges  as  long  7  is  greater  than  7' . 
In  the  integral  expression  (3.24)  for  Snj(7) ,  r' is  restricted  to  0<r'<a,  whereas  r  is  always 
greater  than  a.  In  the  integral  we  will  consider  only  the  first  two  terms  in  the  expansion. 


1  \(,  .  r’)  (r  +  r'cos(y)) 


Likewise  the  exponential  term  in  (3.24)  can  be  written, 


(3.28) 


\r-r\~- 


■  ~  r(l - cos(j'))  =  r  -r'cos(y). 


1+ — cos(y) 
r 


(3.29) 


Expressing  7  and  r  in  cartesian  component  form,  an  expression  for  cosine  (y)  can  be 
developed  in  terms  of  y,  r,  y ' ,  and  r '  . 


r  »r 


=  rrcosj'  =  |(rsin^)j'+  (rcos0)k'j*((/-sin0')j  +  (r'sin0)ij  =  |>’j,  +  zk'j*^>'j  +  xij 


COS 


(y)  =  4 


LI 

r'  r 


Substituting  back  into  the  equation  3.24  for  S„, 


(3.30) 
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r 


where  a  =  k„y/r,  y\x)  =  (a2  -x'2)112  (Figure  3.4),  and  “a”  is  taken  to  be  the  propeller 
radius  (2  m).  Now  all  that  remains  is  to  evaluate  each  of  the  double  integrals  Ii  and  I2. 


Performing  the  first  integration  over  y'  on  integral  I, , 
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/•a 


h  = 


•y”(x) 


exp(-i-a  -y')  dy'dx' 


-a  */-y"(x) 


ra 


h  = 


(exp(-i  a  y")  -  exp(i  a  y"))  — —  dx' 


-i-a 


ra 


h  = 


21 dstlia 

a 


-a 


Evaluating  the  I2term  in  the  same  manner, 


(3.31) 


Oy"(x’) 

y'-exp(-i-a  y')  dy'dx' 

_  x/” 


I2  = 


a  ry"(x') 
a  *'-y"(x') 

(i-y"(x)q  +  l) 


-exp(-i-y"(x)-a)  +  exp(i-y”(x)  a)  ^‘  y  — — 

a  a 


dx1 


-a 


I2  =  — •  f  (y"(x)  a  cos(y"(x)-a)  -  sin(y"(x)-a))  dx' 
a  •'-a 


(3.32) 


The  remaining  integral  over  x'  cannot  be  done  analytically  for  either  Ii  or  I2,  so  a 
numerical  method  was  used. 


In  the  same  manner  as  for  Model  A,  MathCAD  was  utilized  to  manipulate  the  arrays  and 
generate  values  for  the  acoustic  intensity  for  each  observation  point.  Since  the 
experimental  data  does  not  indicate  higher  order  harmonics,  only  the  contributions  of  the 
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fundamental  frequency  (BPF)  of  68  Hz  was  considered.  Plots  of  the  model  results  at  the 
three  phase  angle  sets  defined  in  the  Experimental  section  are  shown  in  Figures  3.5-  3.7. 
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Figure  3.5.  Predicted  Fuselage  Acoustic  Intensity  -  Model  B  (Radiating  Disk)  •  Phase  Angle  (180°, 
180°,  180°,  180°) 
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Figure  3.7.  Predicted  Fuselage  Acoustic  Intensity  -  Model  B  (Radiating  Disk)  -  Phase  Angle  (180°, 
180°,  150°,  120°) 


34 


Model  C  -  Uniformly  Rotating  Paddle 


Instead  of  a  pulsating  disk  as  in  Figure  3.4,  envision  the  source  as  a  single  propeller  of 
radius  “a”  confined  to  rotate  about  the  hub  at  a  frequency  (0  equal  to  the  BPF.  From  the 
acoustic  wave  equation  (3.1),  we  identify  the  source  of  the  radiation  as  a  body  force 
density  F(r,t )  which  actually  arises  from  a  paddle-air  contact  force  at  the  paddle  surface. 
For  simplicity,  we  consider  the  contact  force  as  acting  only  along  the  blade  of  the  paddle 
(assumed  infinitesimally  thin),  with  a  magnitude  at  each  point  on  the  blade  proportional 
to  the  blade  velocity,  i.e. 

^Contact  Per  un'1  length  of  the  blade  =  Kv(r)  =  k(cbxr)  =  kcoplj),  0<p<a, 
where  K  is  the  force  coefficient  per  unit  length  and  the  polar  coordinates  p  and  <J>(=eot-8) 
locate  the  point  on  the  blade.  Here  5  is  the  phase  angle  of  the  rotating  paddle.  The  total 
body  force  is  therefore, 

a 

FTo,ai  =  ^(op4>dp  (3.33) 

0 

By  definition,  FTo!al  can  be  written  as  a  volume  integral  over  the  body  force  density 
F(r,t) ,  so  in  cylindrical  coordinates  (p,<t>,z)  we  have, 

«  «o  2n 

F Total  ~  |  J  J pF(r,t)pdpd<fidz  (3.34) 

o-~o 

Equating  the  two  expressions  (3.33)  and  (3.34)  shows  that  F{r,t)  must  have  the 
representation, 

F{r,t)  =  K(v0d{z)6{<fi-wi+6)  it{OSp<a) 

F(r,t)  =  0  if  (p>a) 
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The  delta  function  over  z  restricts  the  body  force  density  to  the  x-y  plane.  The  other  delta 
function  ensures  periodicity  and  eliminates  harmonic  contributions.  Since  F  has  only  an 
azimuthal  component,  the  source  term  for  the  inhomogeneous  acoustic  wave  equation 
(3.2)  is 

VDF(r,f)  =  — — {K0)S(z)8{<p  -cot +  5)) 

P  d<t> 

VCF(r,f)  = — KCt)S(z)S'(<p  -ax +6)  (3.36) 

P 


Substituting  our  source  term  into  equation  3.13, 
"  1  *(0 


A/»(r,/)  =  ttwf  \—,rz^4iS\<p'-cot  +  S)S 
J±P  \r-r] 


t'-t  + 


Integrating  over  t '  is  straightforward. 


AP(r,f)  =  m  f '~f~riS\0'-oxm  +#)dr\ 

Ip  \r-r\ 


(3.37) 


where, 


^_rTh 


=t~ 


V  J 


Integration  over  z'  gives 

a  2  n 

A/>(r.O  =  M»|  J,  .8X<j>'-ax„,  +S)d<t>'dp\  (3.38) 

0  o  *r  r  I 

Before  we  can  proceed  further,  we  express  the  quantity  |r  -?  j  in  the  integrand  in 


cylindrical  coordinates, 


a  7n 


AP(r,o  -  mv  J  J 


0  0  (r2  +  p'2-2p'(xcos(0)  +  ysinty))} 


SX<f>'-cotret  -vS)d<t>'dp\ 
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using  the  component  expansions  of  r  and  r ' , 

7  =  p'  cos((p')x + p'$\n{(p')y 
r  -  xx  +  yy  +  zz* 

We  are  ready  to  integrate  over  <|> '  . 


Integration  involving  a  delta  function  derivative  is  simply  the  negative  of  the  derivative 
of  the  function  evaluated  at  the  argument  of  the  delta  function.  Thus, 


^’(xsin(6»fe,  -S)-ycos(cotre,  -8)) 
r2  +  p’2-2p '(xcos(cotrel  -S)  +  ysin(oxret  - 


dp\ 


If  tret  is  independent  of  p ' ,  then  y  is  also  independent  of  p '  and  the  quantity  in 
parentheses  in  the  numerator  is  a  constant  in  the  integration.  The  denominator  is  merely 
the  distance  between  the  source  and  the  observation  point;  therefore  we  may  rewrite  the 
integral  as, 

AP(?,,)  =  cJ— - — £ - -5 ndp\  (3.39) 

o[r2  +  /J,2-2/j’cos(j')| 

where, 

x  cos(oxret  -  8)  +  y-sin(oXret  -  8) 
cos(v)  - - i - - - - - - 


and 

C  =  K‘(o^x  -cos (ox  fet  -  8j  -  y-sin^ox  ret  -  8jj 
Making  the  substitution  u  =  p'-rcos(y)  gives 

a-rcos()') 

wc  J  M  +  rCOS(y)  3  du 

-rcos()')  (w2+r2sin2(y))2 
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A P(r,t)  =  C[/3  +  rcos(y)/4]. 


where, 


fl-rcos(j') 

^  I 


"  - du= -l 


J  l  r  I  ’ 

rcosir)  ^u2  +  r2sin2(y))2  (r2  +  a2 -2arcos(y)j2 


and 


fl-rcos(>') 


/<=  J 


1  -A-  1 


rcosir)  (M2  +  r2  sin2(x)) 


f~~  r2  sin2  (j^) 


.  .  o-rcos(y) 

cos(y)  + - i - j- 


(r2  -ha2  -2arcos(y)j2  J 


Therefore, 


A P(r,t)  =  C 


r 


[r2  +  a2  “2arcos(y)p 


cos2(y) 

rsin2(^) 


cos(^)  (fl-rcos(p)) 
rsin2(y)  (r2  +a2  -2arcos(y)f 


AP(r,t)  =  Ccse2(y) 


r 


1 — cos(y) 

_ r _ 

I 

(r2  +a2  -2arcos(y))2 


(3.40) 


The  gives  the  acoustic  signal  from  the  engine  at  each  observation  point  at  time  t.  Each 
engine  will  have  a  different  location  and  a  different  phase  angle. 


Using  the  index  “f’  to  denote  each  engine,  the  sound  level  measured  at  the  BPF  at  each 
field  point  is  proportional  to, 


(3.41) 


In  the  same  manner  as  Models  A  and  B,  MathCAD  was  utilized  to  manipulate  the  arrays 
and  generate  values  for  the  acoustic  intensity  for  each  observation  point.  Plots  of  the 
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model  results  at  the  phase  angle  sets  defined  in  the  experimental  section  are  shown  in 
Figures  3.8-3.10. 
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Figure  3.8..  Predicted  Fuselage  Acoustic  Intensity  -  Model  C  -  Phase  Angle  (180°,  180°,  180°,  180°). 
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Figure  3.9.  Predicted  Fuselage  Acoustic  Intensity  -  Model  C  -  Phase  Angle  (120°,  150°,  180°,  150°) 


Figure  3.1011.  Predicted  Fuselage  Acoustic  Intensity  -  Model  C  -  Phase  Angle  (120°,  150°,  180°,  180°) 
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Section  4  -  DISCUSSION  AND  CONCLUSIONS 


Model  A  appears  to  exhibit  a  couple  of  characteristics  of  the  actual  data  set.  Phase  angle 
set  180°,  180°,  180°,  180°  (case  1)  produced  an  intensity  region  in  the  plane  of  the 
propeller.  This  result  correlates  well  with  the  actual  data  except  that  the  observed  data  has 
the  band  shifted  slightly  aft;  a  phenomenon  that  seems  largely  attributed  to  the  speed  of 
the  aircraft.  A  relative  observation  point  was  determined  based  upon  the  amount  of  time 
required  for  the  wave  generated  from  each  source  to  propagate  to  the  y-coordinate  of 
each  observation  point.  These  changes  in  the  path  distance  corresponded  to  a  wider 
intensity  band  that  included  the  region  indicated  by  the  observed  data.  By  contrast, 
comparison  of  the  120°,  150°,  180°,  180°  (case2)  and  120°,  150°,  180°,  150°  (case3) 
phase  sets  did  not  compare  well  with  the  observed  data  sets.  Variations  in  the  relative 
propeller  phase  in  Model  A  did  not  result  in  the  unusual  lower  intensity  zone  on  the  RHS 
of  the  aircraft  as  in  Figure  2.7  nor  the  large  cool  regions  identified  in  Figure  2.8,  only 
variations  in  the  width  of  the  propeller  plane  intensity.  Model  B  did  not  correlate  with 
the  actual  data  to  any  greater  degree  than  Model  A. 

Model  B  generated  plots  that  exhibited  a  considerable  amount  of  detail.  The  distributions 
generated  for  case  1  match  the  data  set  fairly  well  with  the  exception  that  the  regions  of 
high  acoustic  intensity  and  low  acoustic  intensity  were  reversed.  The  actual  data  set  had 
higher  levels  of  noise  in  the  BPF  near  the  perimeter  of  the  fuselage,  whereas  Model  B 
predicted  lower  noise  levels.  Analysis  of  the  other  two  cases  did  not  match  up  nearly  as 
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well.  Modeled  results  for  case  2  and  case  3  produced  distributions  that  were  fairly 
similar  to  each  other;  an  area  of  low  intensity  along  the  RHS  of  the  graph.  Model  B  did 
not  recreate  the  shift  in  the  low  intensity  region  from  the  RHS  to  the  LHS  of  the  aircraft, 
as  the  phase  on  engine  1  is  incremented  by  30  degrees  from  150  degrees  to  180  degrees 
nor  did  it  predict  rather  large  region  of  low  intensity  exhibited  by  case  3.  Variations  in 
the  relative  phase  between  the  engines  produced  noise  zone  shifts  in  the  z  direction  , 
down  the  fuselage. 

Model  C  performed  equally  as  well  as  Model  A  in  producing  a  region  of  high  intensity  in 
the  propeller  plane  when  the  engines  were  in  phase,  which  corresponds  fairly  well  to  the 
observed  data.  However,  as  in  the  other  two  models,  variations  in  propeller  phase  did 
not  recreate  the  shift  in  the  noise  field.  In  fact,  as  the  engine  #1  cycled  from  case  2  to 
case  3  the  predicted  noise  levels  actually  increased  in  the  propeller  plane.  Phase  changes 
contributed  to  lateral  shifts  in  the  distributions  in  the  propeller  plane. 

This  thesis  developed  a  framework  for  determining  the  acoustic  intensity  at  an 
observation  point  based  upon  Green’s  Function  techniques.  Three  models  were  initially 
considered,  a  monochromatic  point  source,  an  extended  pulsating  disk,  and  a  uniformly 
rotating  paddle.  Although  all  three  exhibited  some  characteristics  of  the  actual  data  set, 
they  ultimately  failed  to  adequately  predict  the  shifting  in  the  noise  field  associated  with 
controlled  variations  of  the  relative  propeller  phase.  The  actual  source  function  may  be  a 
combination  of  Model  B  and  Model  C.  i.e..  some  sort  of  extended  line  source,  since 
Model  B  produced  variations  in  the  z  direction  and  Model  C  the  y  direction.  Future  work 


44 


should  center  on  developing  more  appropriate  functions  to  represent  the  source  and 
investigate  application  of  these  techniques  to  the  far-field  noise. 
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